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Abstract 

We suggest a new approach for calculating heavy-to-light form factors. The method 
is based on light cone sum rules (LCSR) and covers the whole kinematical range of 
momentum transfer. The derivation of the new sum rule uses a suitable combination 
of double and single dispersion integrals. As an example we give numerical results 
for the form factor for the i? — > vr transition. 
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1 Introduction 



The accurate study of B meson decays is a main source of information for understanding 
CP violation. We expect that the upcoming experiments will measure a variety of B 
decay properties with good precision [|l|. In order to over-constrain the unitary triangle 
one pursues not only a measurement of the various angles, but tries as well to determine 
the length of the sides of the triangle. One important quantity is the CKM matrix element 
\Vub\, which is proportional to the length of one side of the unitary triangle. \Vub\ can 
be obtained from semileptonic B decays. There are two complementary strategies for a 
determination of \ Vub\, relying either on exclusive or inclusive measurements. The inclusive 
decay B Xlu can be calculated using the heavy-quark expansion |[. However one 
of the main obstacles here is that experimentally it is mandatory to impose restrictive 
cuts to suppress the background from B — >• Xclv decays. 

On the other hand an exclusive measurement through the channels B ^ nlu or B ^ 
plv is easier from an experimental point of view. It requires however a reliable theoretical 
calculation of the heavy-to-light transition form factors. These form factors parameterize 
the relevant hadronic matrix elements and are non-perturbative quantities. They can be 
calculated using lattice methods or QCD sum rules on the light-cone. In this paper 

we focus on the sum rule approach. 

The QCD light-cone sum rule method (LCSR) has been suggested in [10|, [Tl[| and 



is a combination of the operator product expansion (OPE) on the light-cone [|l^, 0, |14 



with QCD sum rule techniques [|T5|. For a review of the method and results we refer the 



reader to 17 



Within the sum rule approach results for the heavy-to-light form factors corresponding 
to the B^Ti, B^K,D^7i or D—>-K transitions are easily obtained from one master 
formula ||T^-|2^. Changing the light meson from a pseudo-scalar meson to a vector meson 
(e.g. p, uj, K* or 0) the form factors of semi-leptonic or rare radiative decays of a 5-meson 
into vector mesons are obtained |]2S|, ^ ^ . 



To illustrate our method we focus on the form factors relevant to the B n transition, 
defined through 

{7r{q)\u^,b\Bip + q)) = 2/+(/)g^ +(/+(/) + /-(/)) (1) 

where p + q, q and p denote the B and vr four-momenta and the momentum transfer, 
respectively, and are the two independent form factors. If one neglects lepton masses 
only the form factor /"*" is relevant. The standard sum rule technique starts from the 
correlation function of two heavy-light currents 

F^ip,q) = 1 1 dxe'P-'' <n{q)\T{u{xh^b{x),mbb{0)n5d{0)}\0> . (2) 

The correlation function is then expanded near the light-cone ~ 0. Each power of in 
the light-cone expansion leads to an additional power of ml — {p+uqY in the denominator. 
Here u denotes the momentum fraction carried by one of the quarks inside the pion. The 
light-cone expansion is justified, provided that {p + qY and p^ are sufficiently smaller than 
ml. In particular this implies that p^ < m^ — 2mbX, where x is a scale of order Aqcd- 
This limits the application of the sum rule to the kinematic region to 

< < - 2mbX- (3) 
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However, in order to extract \Vub\ one needs a prediction for the form factor /"•" in the 
whole kinematical range < < {rriB — fnT^Y, which is much wider than (j^). Several 
attempts have been made in the past to overcome this problem. One of the first solutions 
23| consisted of extrapolating the form factor from the region < < — 2mhX to 



the region ml — 2mi,x < < {ms ~ Tn^^Y by assuming the functional form 
Another analysis used the fact that close to the point p"^ = m\, the form factor can be 



described by a simple pole model (see for example [28, 29, [S^l), assuming vector meson 
dominance: 

f+f„2N _ fB*gB*Bn , . 

^ ' 2mB* {1 - pyml,) ^ ' 

The coupling gB*Bn can be calculated by QCD sum rules techniques, starting from the 
same correlation function eq. (0) and using a double dispersion relation. This additional 
information was used in by fitting a function of the form of eq. (|) to the low p^- 
behavior obtained from the standard light-cone sum rule and to a single pole model 
description for high p"^. A third approach was used in the analysis of ||22|: The Becirevic- 
Kaidalov parameterization |^T| 

^^^^'^ " '''' [l-pyml, ~ l-pV(7m|.)) 

was assumed as functional form and the parameters of this ansatz were determined from 
QCD sum rules, using the light-cone sum rule for < < ml — 2mbX as well as 
the information on the coupling Qb'b-k- In the parameterization eq. (|^) the first term 
corresponds to the vector meson pole, whereas the second term represents an effective 
contribution from all higher resonances. 

All of the methods mentioned above involve ad hoc assumptions, which are difficult 
to justify from first principles. In addition we cannot easily quantify the error associated 
with these additional assumptions. 

In this letter we suggest a new method for calculating heavy-to-light form factors which 
yields a prediction in the whole kinematical range of momentum transfer and overcomes 
the problem outlined above. We start from the same correlation function eq. (0). The 
derivation of the new sum rule is based on a combination of double and single dispersion 
integrals. The particular combination of double and single dispersion integrals ensures 
that our sum rule is valid over the whole kinematical region. As input data we need the 
value of /^(O), together with the first / derivatives /^'•'•'(p^) at = (with / an integer), 
which can be obtained (numerically) from the standard sum rule for f~^{p'^). In addition 
we need the value of the coupling gB*B-ni which can be obtained from the sum rule for 
the coupling. We derive the new sum rule to the same accuracy to which the two other 
sum rules are known I^O], |2^ : to next-to- leading order in twist 2 and to leading order 



in twist 3 and 4. We present numerical results for the form factor f^{p'^) for the B ^ n 
transition. 

This paper is organized as follows: In the section 2 we introduce the new sum rule. 
In section 3 we give for the new sum rule the QCD corrections to twist 2 relevant for 
the B ^ 71 transition. Numerical results are given in section 4. Section 5 contains our 
conclusions. 



3 



2 The formalism 

We write the correlation function eq. (|^) in terms of invariant amplitudes {p + q)^) 

and F{p^, {p + g)^) 

F^{p,q) = F{p\{p + qf)q^ + F{p\{p + qf)p^, (7) 
and focus on (p + qY). We denote by 

p(si,S2) = -i^ Im^i Im^a F(si, S2) (8) 

the imaginary part of F with respect to S2 and the double imaginary part of F, respec- 
tively. Furthermore we define by 

B,2 = hm _L^(_p2\"My (g) 

p2^-oo,n^oo,-p2/n=M2 [n — 1)! ^ ^ \dp J 

the Borel operator with respect to the variable p^. Hadronic representations for the 
spectral densities o"(p^, S2) and p(si, S2) are written by singling out the ground state: 

a^-^\p\s2) = 2mlfBf^{p')6{s,-ml) + a^^'^{p\s2)Q{s,-s,), 
p'""^'"(si, S2) = m%mB*fBfB*gB*B7TSisi - m|,)(5(s2 - rn%) 

+P^''*^(si, S2) (1 - e(si, S2)) e(si - m2)e(s2 - m^). (lo) 

Here 6(si,S2) defines the duality interval in the (si, S2)-plane for the ground state. A 
convenient choice is given by the square 

e(Si,S2) = 0(So-Si)0(so-S2) (11) 

The standard sum rule for the form factor /^(p^) is obtained from writing a single disper- 
sion relation for -F(p^, (p+q)"^) in the (p+g)^-channel, inserting the hadronic representation 
and Borelizing in {p + q)"^'- 

The Borel operator ensures that any subtraction terms which might appear will vanish 
after Borelization. One proceeds to replace a^"''^ by a^^^ and equate the r.h.s of eq. (p!2| ) 
to the QCD calculation of B(p+q)2F(jP' , {p + g)^). This yields the standard sum rule for 

r(/) = ^^l^''''''(p"^'-y~'^ds, (13) 
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In a similar way the standard light-cone sum rule for the coupling Qe'e-k is obtained from 
a double dispersion relation: 



m%mB'fBfB*gB*B-n 



+ / dsids2-, 



p^)(s2 - (p + qy) 



2^ ' 



(14) 



where S denotes the union of Si > Sq, S2 > ml with si > ml, S2 > Sq. Again the presence 
of the Borel operators ensures that any subtraction will give a vanishing contribution. 
The LCSR for the coupling gB*B-K reads: 



gB*BiT 



m%mB'fBfB' 



sp So 

dsi / c?S2p'^^^(si,S2)e" 



•B* ' "-"'^ ' B' 



(15) 



9 9 

mf mf 

b h 



Here we took the two Borel parameters to be equal. 



To derive our new sum rule we suggest here to use a dispersion relation for cr(p^, S2)/(p^)' 
in the p^-channel (with I being an integer) 



ct(p^ S2) 



1 / 2\' d'- ^ cr(si,S2) 



P 



(/- 1)! V" ^ ds[-' si-p 



+ / ds^J^'-^ (16) 
Si Si — p'^ 



and to replace aijP', S2) in eq. (|^) by the r.h.s of eq. (|^). By choosing / high enough the 
dispersion relation eq. (|TBp will be convergent. We obtain 

%+,)2F(p2,(p + g)2) = 

^2m|/s/+(p2) r ^^^^^^(p5^_ P{si,S2) 

J ^ ^ s{ {Si 



m| - (p + g)" 



v 



si>m^,S2>so 



p'^){s2 - {p + qY) 



ds2 ^^2\i d'' ^ o-{si, S2) 



(/-I)! J S2-{p + q)^''^ ' ds[-' si-p 



S2>S0 



(17) 



si=0/ 



Furthermore we write down a double dispersion relation for F ijp' , {p + q)"^) / ijP'Y : 
B^,+,)2F{p\{p + qf) = 

{p^Y mj^mB*fBfB*gB*Bn 



Bo, 



{m%,y (p2 - m\,){{p + g)2 - m|) 



dS2 ^^2\l ^ Or(si, S2] 



(p 

+ I dsids2 — ; 



2\l 



P(S1,S2) 



4 (si -p2)(s2 - (p + g)2) 



(/_!)! J s2-{p + qT'^ ' ds'f' si-p 



S2>m 



si=0 



Again, by choosing / high enough the dispersion integral will be convergent in the Si- 
channel. The Borel operator ensures that any subtraction terms in the S2-channel will 
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vanish. Now equating the r.h.s of eq. ([T7|) with the r.h.s of eq. ([T8|) we obtain the sum 
rule 



n+( 2^ ^ 1 {P'^y fB*9B*Bn 1 / 2y ^' ^ f^i^l 



si=0 



where the region S' is defined by si > sq and ml < S2 < sq- This sum rule is valid in the 
whole kinematical range of p^. As input data we need the first (/ — 1) terms of the Taylor 
expansion of f^{p^) around = 0. These parameters can be obtained numerically from 
the standard sum rule eq. ([T^). We further need the residuum at the pole = m^*, which 
can be obtained from the sum rule eq. (pTSf). The new sum rule agrees by construction 
with the standard sum rule eq. (|T3) in a Taylor expansion around = up to the first 



(/ — 1) terms. Furthermore the residuum at p"^ = m\, agrees with the coupling sum rule 
eq. (^. 

We remark that the parameter / plays a similar role as the Borel parameter M ^ : There 
is a lower limit on I since the dispersion relations eq. ( p!6[ ) and eq. (|18|) have to converge. 
Going to higher values for I will improve the convergence of the dispersion relations and 
suppress higher resonances in the i?*-channel. But there is also an upper limit on I: The 
higher the value of Z, the more derivatives of f~^{p'^) at = we have to know. At 
some point we start probing the region p^ > ml — 2xmb, at which the standard sum rule 



eq. (|T3|) might break down. By using the sum rule with various values of /, say I = 1, 2, 3 
and by looking at the variation of the results, we can get an estimate of the uncertainty 
of our method. 

In the case / = the second term is absent. For I = the first term corresponds to 
the pole model. As we will show explicitly below, the last term vanishes for the leading 
order twist 2, 3 and 4 contributions. The first non-vanishing contribution comes from the 
ttg-corrections to the twist 2 contribution. This might explain the empirical fact that for 
some heavy-to-light transitions (like for D ^ tt) the simple pole model approximates the 
form factor reasonably well. 



3 The additional term 

We now consider the last term in eq. 




fLAp-) = f^s,,s,_ ( t] P(£^e-^. (20) 



which has to be evaluated. The sum rules eq. (|T3D and eq. (|T5D are known in twist 2 to 
NLO accuracy, and in twist 3 and 4 to LO. Aiming at the same accuracy for eq. ( PO]) we 
find that the LO contributions of twist two, three and four vanish. This is due to the fact 
that the spectral density p{si, S2) at leading order is localized along the diagonal si = S2, 
whereas the integration region E' lies beyond the diagonal. At next-to-leading order the 
situation is different. The radiative corrections give contributions which are smeared over 
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the whole Si — S2 plane, overlapping with the region S'. For the twist 2 NLO contribution 
we find: 



where 



— / dr- 



9i[r) 



92[r) 



r - If 



2 

— - I dr 



exp 



mt — m 



B 



M2 



dr- 



9i{r) 
(r-1) 



9_m 

r — 1 



+ 2 



dr 



g2{r) ln(r) 
(r - 1)3 



~ I dr 



r — 1 



(l+r)^o 



ds r e n 



+ In 



1 + r 



1 



,1 + r + s 

(l+r)zo 

ds r e !+'■ 



1 + r + rs/ 1 + r + rs — a — ar 

1 + r + rs\ \ 
1 + r + s // ' 

1+r \' 1 



1 + r + rs/ 1 + r + rs — a — ar 



:2i) 



:22) 



Here we used the dimensionless variables 



m 



2 ' 



So - 



mt 



(23) 



Eq. (^) is written in such a form that singularities at r 1, which are present in 
individual contributions, cancel explicitly in the combinations inside the square brackets. 
The values of gi and g2 and their derivatives at r = 1 are 



9iil) 
9'lil) 

92{l) 



0, 9'iil) 



1 



1 



(1 + 2:0)' I + ZqI + Zo - a' 



—e 



-bzo 



1 



1 



■ (j,Z()a — 2^00 — bzQa 



(1 + ZoY (1 + zoy {l + zo- ay 
—bz^a — a + bzo + bz^ + 2bzl + z^ — Izq — Izq + 1 + 2zq 







^2(1) = -e 



(1 + zoY 1 + zo - a' 
The integrations over r and s can be performed numerically. 



(24) 



4 Numerical results 

We perform a numerical evaluation of the new sum rule eq. (0) for the values / = 
0, 1, 2 and 3. We need therefore the values of /"*" and its first two derivatives at p"^ = 0. 
These numbers can be obtained from the sum rule eq. (0). The derivatives are obtained 
numerically according to 
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/t(o) 


(/^.)'(0) 


(/^.)"(0) 


0.28 


0.014 GeV"" 


0.0014 GeV"* 



Table 1: The form factor f^^ and its first derivatives at = GeV. 



The sum rules eq. (1T3D and eq. ([T5| ) depend on various input parameters, where each 
parameter is only known within a certain range. A complete error analysis, where each 
parameter was varied within an interval, was carried out in In this letter we do not 



repeat such a complete error analysis. We only study the dependence on I and we fix 
here these additional parameters to certain values. As numerical input parameters we use 
rrib = 4.7 GeV, Sq = 35 GeV^ and a{mz) = 0.118. We note that the value of the pole 
mass is compatible with the recently updated MS mass rh{rh) = 4.2 ±0.1 (for a review of 
recent results see ) . The same values of mf, and Sq have been used in the two-point sum 



rule for estimating the decay constant /b- We obtained /b = 183 MeV. The coefficients 
of the leading twist pion distribution amplitude have recently been updated in |3^ , 



the coefficients of the twist three and four amplitudes can be found in |jTT], |T8|, |1^, |35 . 
Table (|I|) shows the results for the form factor and its first two derivatives at = 
calculated from the sum rule eq. (O). We further need the 3*871 coupling as input data. 



This value is obtained from the sum rule eq. (0) as fB*9B*BTr = 4.4 GeV [|TU|, ^ 



With these input values we evaluate the sum rule eq. ([19D for / = 0,1,2 and 3. For the 
Borel parameter we use M2 = 10 GeVl Our results are shown in fig. (|l|). First of all we 
note a remarkable stability of the numerical results with respect to changing the number 
of subtractions /. The results for / = 1, 2 and 3 are almost identical. Secondly, the result 
for / = (which is the pole model plus an ag-correction term) differs at low momentum, 
but approaches for high momentum the results with subtractions. This corresponds to 
the known fact, that for the B tt transition the pole model does not describe the form 
factor accurately at low momentum. This is also shown in fig. (|^). We also note that 
the standard sum rule for /"*" eq. (|T3) will differ significantly from our results at high 



momentum. Our results for I = 1, 2, 3 agree well with the parameterization given in |22 



/+,(0) = 0.28 ± 0.05, as. = 0.32 ±°:^^ . (27) 
To see the difference we take / = 2 as our main result and plot in fig. (0) the deviation 



9ip ')-fU p'' 

fUp' 



^ = — 7^-7;;^^ — 



where we take for g the results fiL-^^, f,t^ and the values according to eq. (^) with 
/^^(O) = 0.28 and = 0.32. The deviations are small. 

In fig. (^) we show the sum rule results for the B ^ n form factor, eq.(|l9D with 
/ = 2 (solid line) and the pole model prediction (dotted line) in comparison to lattice 
results. The lattice results come from FNAL ||^ (full circles), UKQCD ||^ (triangles), 
APE [§ (full square), JLQCD § (open circles), and ELC § (semi-full circle). Taking 



into account the uncertainty of our result, which is estimated to be roughly 15 — 20% |22 



we observe a satisfactory agreement with most lattice results. However, the results from 
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the JLQCD collaboration lie systematically above our values. In addition, this group 
quotes rather small uncertainties. In general, further improvements in the accuracy of 
lattice calculations are welcome. 



5 Conclusions 

In this letter we have shown a method how to obtain the form factor for heavy-to-light 
transitions in the whole range of momentum transfer. Our method extends the QCD sum 
rule approach and uses a combination of single and double dispersion relations. It involves 
an additional (integer) parameter I, corresponding to the number of subtractions in one 
channel. We have derived the corresponding sum rule for the B —>■ n transition to twist 
four accuracy and including radiative corrections to the twist-2 contribution. As input 
data we need the first (Z — 1) terms of the Taylor expansion of f~^{p'^) around = as 
well as the residuum at the pole = m^, , which can reliable be obtained from standard 
sum rules. The new sum rule involves an additional term, which vanishes in leading order 
for the twist two, three and four contributions. We have calculated the non-vanishing 
next-to-leading order twist two contribution to this term. We have shown that variation 
of the parameter / introduces only small numerical changes in the final result for the form 
factor, which are negligible against other uncertainties. 
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Figure 1: The results of the new sum rule for the form factor fB^niP^) fo^' different 
numbers I. The dotted line corresponds to Z = 0, the solid hues correspond to / = 1, 2, 3. 



0.1 




-0.05 ' ' ' ' ' ^ 

5 10 15 20 25 

Figure 2: The relative deviation of the the numerical results with / = 1 from / = 2 (solid 
line), of / = 3 from / = 2 (dashed line) and of the result quoted in from / = 2 (dotted 
line). 
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Figure 3: The sum rule results for the B ^ n form factor, eq. ([l9|) with / = 2 (solid 
line) and the pole model result (dotted line) in comparison to lattice results. The lattice 
results come from FNAL § (full circles), UKQCD § (triangles), APE § (full square), 
JLQCD (open circles), and ELC |^ (semi-full circle). 
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